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This article studies , by computer simulations, the performance of deep-space 
telemetry signals that employ the pulse code mod ulation /phase modulation 
(PCM/PM) technique, using nonreturn-to-zero data, under the separate and com- 
bined effects of unbalanced data, data asymmetry, and a band-limited channel. The 
study is based on measuring the symbol error rate performance and comparing the 
results to the theoretical results presented in previous articles. Only the effects 
of imperfect carrier tracking due to an imperfect data stream are considered. The 
presence of an imperfect data stream (unbalanced and/or asymmetric) produces un- 
desirable spectral components at the carrier frequency, creating an imperfect carrier 
reference that will degrade the performance of the telemetry system. Further dis- 
turbance to the carrier reference is caused by the intersymbol interference created 
by the band-limited channel. 


I. Introduction 

There is considerable interest among international space agencies in searching for a bandwidth-efficient 
modulation scheme that can be used for future space missions without major modifications to their ground 
stations |1-4J. The Consultative Committee for Space Data Systems (CCSDS) has undertaken the task of 
investigating a modulation scheme that offers both of these features (bandwidth efficiency and no major 
hardware modifications to the current systems). 

Currently, the space telemetry systems employ residual carrier modulation with subcarriers that are 
used to separate the data from the RF residual carrier. This was necessary to avoid interference because 
most of the data power fell within the bandwidth of the carrier phase-locked loop (PLL), as shown 
in Fig. 1(a). The CCSDS has recommended that square- wave and sine-wave subcarriers be used for 
the deep-space and near-Earth missions, respectively [5]. This modulation scheme is called pulse code 
modulation/phase-shift keying/phase modulation (PCM/PSK/PM), and it was developed at a time when 
weak signals and low data rates dominated [6], With the development of technology and the evolvement 
o t e Deep Space Network (DSN), a significant increase in the signal power can result in higher data 
mtes. sing subcarriers in this case causes the occupied bandwidth to increase significantly. This is 
prohibitive because the space telemetry systems often operate under imposed bandwidth constraints. A 
natural solution is to eliminate the subcarrier and modulate the nonreturn-to-zero (NRZ) data directly 
on the RF carrier. This modulation scheme is referred to as PCM/PM/NRZ, and not only does it require 
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Fiq. 1. PCM/P M/NRZ modulation: (a) low data rate: high 
ratio of loop bandwidth to data rate and (b) high data rate: 
low ratio of loop bandwidth to data rat©. 


minimum hardware modifications to the current systems, but it also achieves the bandwidth efficiency 
[4 71. In this modulation technique, the part of the data spectrum that falls within the narrow carrier 
loop bandwidth seems flat and appears as white noise, as shown in Fig. 1( ), an smce e ra 10 o 
bandwidth to data rate is very small, the carrier tracking performance degradation due to this white 

noise component is negligible. 

Recently, Nguyen has investigated and analyzed the behavior of PCM/PM receivers in nonideal chan- 
nels [1,2]. The imbalance between +l’s and -l’s and/or data asymmetry in the data stream produce 
undesirable spectral components that degrade the performance of the system. Further degradation s 
caused by the intersymbol interference (ISI) created by the band-limited channel. This article verifies 
by computer simulations, the theoretical results presented in [1,2] for the NRZ data st rearm The Signa 
Processing Worksystem (SPW) was used for implementing and simu atmg the system. The se g 
effects of unbalanced data, data asymmetry, and band-limited channel on the symbol error rate (SE R) 
performance of PCM/PM/NRZ receivers were simulated and then compared to the theoretical res 
m^nTed Tn [ll. In reality/however, the receivers operate in the aggregate presence of these three effects, 
Ld the symbol signal-to-noise ratio (SSNR) degradation due to the three effects is not the algebraic 
sum of the SSNR degradation due to each separate effect. The second part of this article P^sents the 
simulation results for the degradation due to the combined effects on the SER performance, and these 
results are compared to the theoretical results presented in [2]. 

The organization of this article is as follows: Section II describes the separate effects on PCM/ p M/NRZ 
receivers of perfect, unbalanced, asymmetric, and band-limited data streams. Section III describes th 
combined effects of these streams on PCM/PM/NRZ receivers. Section IV gives a brief description of 
the PCM/ PM receiver blocks that were used to build the system in the SPW, Section iscusses e 
simulation results and compares them to theory, and, finally, Section VI presents the conclusion. 


II. Separate Effects on PCM/PM/NRZ Receivers 

The deep-space received telemetry signal, in the absence of a subcarrier, is given by [1] 
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s r {t) = V2P {cos (mr) cos (w c t + do) - d(t) sin sin (u; c t 4- do)} + n(t) 


( 1 ) 


where P is the transmitted power, m T is the telemetry modulation index in rad, u) c = 2i r/ c is the angular 
carrier center frequency in rad/s, #0 is the initial carrier phase, n(£) is an additive white Gaussian noise 
(AWGN), and d(t) is the data stream (NRZ) defined as 


00 

d{t) = ^2 d kP(t + kT s ) (2) 

fc=— OO 

where d* = ±1 with transition density p t , p(t) is the baseband pulse, and T s is the symbol period in 
seconds. The first and second terms of Eq. (1) are the residual carrier and data components, respectively. 

The undesired spectral components caused by the imperfect data stream (unbalanced data and/or 
data asymmetry) can degrade the carrier tracking performance. If d denotes the carrier loop estimate 
of do, the phase error due to the thermal noise and interference caused by the imperfect data stream is 
defined as 


@e = @0 — d = O e (noise) + 9 e (data) + 9 e (spike) (3) 

where d e (noise), 9 e (data), and 9 e (spike) are the phase error caused by the noise, data interference, and 
the spike caused by the imperfect data stream, respectively. 

The carrier loop tracks the residual carrier component in Eq. (1) to provide an imperfect reference 
given by 


r(t) = y/2 cos (a ) c t + d) 


( 4 ) 


The average probability of error due to the imperfect carrier tracking is given by 

Pe = J Pe(0 e )P(0e) d6 e 


( 5 ) 


where P e (0 e ) is the conditional probability of error, and P(9 e ) is the probability density function (pdf) of 
the carrier tracking phase error d e . Assuming that this pdf has a Tikhonov distribution that is entirely 
characterized by the mean (assumed 0) and variance a 2 of d e , and when the loop signal-to-noise ratio 
(SNR) is high, P{9 e ) may be approximated as Gaussian distribution, namely, 


m) 


ex P (~ fle/( 2cr2 )) 

[27TCT 2 ] -1 / 2 


, — OO < d e < OO 


( 6 ) 


As mentioned above, this expression was derived assuming the mean of the phase error 6 e to be zero. 
This assumption, however, is not true for an imperfect data stream, as will be shown in the subsequent 
sections. 

The expressions for P e (9 e ) and the carrier tracking phase error variance a 2 have been evaluated in [1,2] 
for all the different cases studied in this article. The final results will be presented here for completeness. 
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A. Perfect Data Stream 

In a perfect purely random data stream, the probability of transmitting a +1 pulse (or probability of 
mark), p, is equal to the probability of transmitting a -1, q, with transition density, p t , given by 


1 

Pt = 2pq = - 


( 7 ) 


where q = 1 — p. 

The carrier term of Eq. (1) generates a residual carrier at f c with power, P c , given by 

P c = P cos 2 (hit) (8) 

Combining the carrier and data terms, the one-sided power spectrum of a PCM/PM/NRZ perfect data 
stream is given by 

S(f) = S c (f) + S D (f) (9) 

where 

Sc(f) = PcS(f) ( 10 ) 


and 


Sp{f) = P sin 2 (m T )S CO nt(f) 
is the data spectrum with power Pp defined as 


( 11 ) 


P D = P sin 2 (mp) 


(12) 


For a perfect NRZ data stream, S cont {f) is defined as the power spectral density (PSD) for an ideal NRZ 
data stream and is given by 


Scont(f) = T s 


j sin 2 (jrfTs) \ 

l («fTs ) 2 J 


(13) 


Figure 2(a) shows the power spectrum of a perfect NRZ data stream (generated using SPW for symbol 
rate R s = 1 /T s = 10 4 kbits/s). 

For a perfect data stream and ideal channel, the conditional probability of error is given by 


Pe{0e) = \ erfc 



(14) 


where E s /Nq denotes the SSNR, that is, 
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Fig. 2. Spectrums of different NRZ data streams: (a) 
balanced data stream (p = 0.5), (b) unbalanced data 
stream (p = 0.4), and (c) asymmetric data stream (£ a 14 
percent). 


E, _ PT a sin 2 (m T ) _ PdT s 

N 0 ~ N 0 ~ N 0 

and erfc (x) is defined as the complementary error function given by 


X 

erfc (x) = 1 - erf (x) = 1 [ exp (~v 2 ) 

V™ J 


Note that for this case, the mean of the phase error 0 e in the steady state is zero. This, however, is not 
true for an imperfect data stream, as will be shown in the subsequent sections. 

For the high data rate case ( B L /R S <<0.1, where B L denotes the one-sided loop bandwidth), the 
variance of the carrier tracking phase error is given as [1] 


o 2 — b tan 2 (m T ) 

Po 
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where 


( 18 ) 


(EJNq) 

Po (B l /R s ) tarn 2 {tut) 


By substituting Eqs. (6) and (14) into Eq. (5) and performing the numerical integration, the curve for the 
probability of error versus SSNR was obtained and is shown in Figs. 3 through 9 for comparison purposes. 
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Fig. 3. Theory and simulation versus SSNR for 
unbalanced data without modification to the phase 
error. 


B. Unbalanced Data Stream 

The imbalance between +l*s and -l’s in the data stream causes an additional corruption to the 
received signal in Eq. (1), generating undesirable spectral components that can potentially degrade the 
performance of the telemetry system. When p is not equal to 0.5 (and, therefore, p t < 0.5), the data 
component will be affected and Eq. (11) now becomes 


So{f) = Psin 2 (m T ) {Sddf) + S cont (f)} 


(19) 


where S dc (f) is the dc (or harmonic) component caused by the imperfect data stream that falls on the 
RF carrier. 

The spectrum of an unbalanced NRZ data stream for p = 0.4, generated using the SPW, is shown in 
Fig. 2(b). For a PCM/PM/NRZ unbalanced data stream, the dc and continuous PSD components are 
found to be [1] 
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Fig. 4. Theory and simulation SER versus SSNR 
for unbalanced data. 



Fig. 5. Theory and simulation versus SSNR for data 
asymmetry. 
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Fig. 6. Theory and simulation SER versus SSNR 
for band-limited channel. 
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Fig. 7. Theory and simulation SER versus SSNR for 
various values of unbalanced data. 
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Fig. 8. Theory and simulation SER versus SSNR for 
various values of data asymmetry. 
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Fig. 9. Theory and simulation SER versus SSNR for 
various values of band-limited channel. 
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Sdc(f) = (1 - 2 p)HU) 


( 20 ) 


Scontif) = 4 T a pq 


( sin 2 (tt/T,) 

l 


} 


with power given by 


P<tc = (1 - 2 p) 2 P sin 2 (m T ) 
Pcont = 4pqP sin 2 (m T ) 


respectively, and where 


Pd — Pdc “b Pcont 


( 21 ) 


( 22 ) 

(23) 


(24) 


Therefore, in addition to the tone generated at f c by the residual carrier component in Eq. (1) with power 
given by Eq. (8), the spectrum of unbalanced PCM/PM/NRZ will include another tone at f c generated by 
the imbalance between +Ts and -l’s with power given by Eq. (22). However, these two tones at f c have 
noncoherent phases, causing the mean of the carrier tracking phase error in the steady state to deviate 
away from zero. This deviation is defined as 9 e (mean) y which is a function of p and the modulation index 
rriT and is given by 


9 e (mean) = - tan 1 {(tan m T )(2p - 1)} (25) 

Note that when p — 0.5, then tan -1 0 = 0, independent of m T , as one would expect. 

Figure 10 shows the 9 e (mean) of balanced and unbalanced data streams as generated by an SPW for 
0o = 0- Note that for the case of a balanced data stream, the mean of the phase error is centered at zero, 
whereas for an unbalanced data stream, the mean is at a negative value which, using the above equation, 
is calculated to be about -0.54 rad (-31 deg) for p = 0.6 and m T = 1.25. 

The conditional probability of error P e {9 e ) is the same as the one given by Eq. (14). Recall, however, 
that Eq. (6) for the pdf of the carrier tracking phase error P(9 e ) was derived assuming the mean of 9 e to 
be zero. Therefore, the simulations will have to compensate for the phase difference (0 O - 9) (Eq. (3)) by 
adding the value of 9 e [mean) (Eq. (25)) to the phase of Eq. (4), (#), which results in a zero-mean phase 
error. In that case, P(9 e ) is given by Eq. (6) with the tracking variance given by 

2 1 CK n , . 1 I 

° = 7 0 + 2 tan (mT)+ 2C (26 > 

where p 0 is defined as before, a is the interference due to the continuous spectrum, and I/C is the 
interference caused by the dc component- to-carrier power ratio given, respectively, by 

a = 47>(1 - P) J \H(2nf)\ 2 df (27) 

-OO 
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Fig. 10. 6e (mean) of balanced and unbalanced data 
streams: (a) balanced data stream (p = 0.5): mean of 
phase error is centered at zero and (b) unbalanced data 
stream (p = 0.6): mean of the phase error deviates away 
from zero. 


I 

c 


— (1 - 2 p) 2 tan 2 (uit) 


(28) 


where H(j2nf) denotes the carrier loop transfer function, and for a second-order PLL is given by 


\H(j2nf)\ 2 


1 + 2 (///„) 

1 + (///n) 


(29) 


where /„ is the loop natural frequency. 

The plot of the SER versus SSNR is shown in Figs. 3 and 4. 

C. Data Asymmetry 

Data asymmetry, due to rising and falling voltage transitions, causes undesirable spectral components 
that degrade the performance of the space telemetry system. The data asymmetry model adopted in this 
article assumes that +1 symbols are elongated by (AT s )/2 (relative to their nominal value of T s seconds) 
when a negative-going data transition occurs, and — 1 symbols are shortened by the same amount when a 
positive-going data transition occurs. Otherwise (when no transitions occur), the symbols maintain their 
nominal T s seconds width. This model is illustrated in Fig. 11 for a purely random NRZ data stream. 
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Fig. 11. Perfect and asymmetric NRZ data 
streams: (a) £ = 0 (perfect data stream) and (b) ^ 

£ = 0 (asymmetric data stream). 

The power spectrum of an asymmetric NRZ random data stream with equiprobable symbols (that is, 
p = p t = 0.5) and symbol rate R s of 10 4 kbits/s is shown in Fig. 2(c). The dc, continuous, and harmonics 
PSD components are given by [1,3] 


Sdc(f) = ^(/) 


(30) 


Scont{f) = Y 


sin 2 {irfT a ) 

. (tt/T *) 2 . 


[3 + 5 cos' (,/r.Ol + S Kjyy -i] [3 cos' (*/r,) + cos' (2nfT,()] 


(*7^) 2 J 


(31) 


1 00 1 1 

Sk(f) = yE ^2 c ( m ' r° 6{f ' mRa) 


(32) 


m= 1 


respectively, where £ denotes data asymmetry and is defined as 


£ = 


A 

2 


and where 


C 



- sin 2 (2m7r£) 
4 


(33) 


(34) 


Hence, the data spectrum can be written as 
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So{f) = P sin 2 (m T ) {S dc (f) + S cont {f) + S h (f)} 


(35) 


which, in addition to the tone at f c caused by the carrier component, generates a spike at f c due to 
the dc component Sdc(f), and a spike at integer multiples of the symbol rate R s due to the harmonics 
component 5^(/). The continuous spectrum S cont (f) is plotted in Fig. 12 for various values of £. Note 
that when f = 0, the above equation reduces to the perfect NRZ random data case given by Eq. (11). 

Similar to the unbalanced data case, the phase of the dc component at f c caused by asymmetry and 
the phase of the carrier tone are noncoherent. The mean of the phase error 9 e for a perfectly balanced 
asymmetric data stream was derived to be 


9 e (mean) = - tan 1 |(tan mr) j (36) 

Note that when £ = 0, Q e (mean) = 0, as expected. Again, the simulations may have to compensate for 
the phase difference (Eq. (3)) to make the mean of 9 e zero at steady state. 

Recall that in order to calculate the average probability of error, the conditional probability of error 
P e (9 e ) and the tracking variance a 2 must be determined. For the data asymmetry model used in this 
article and for a purely random and equiprobable (perfectly balanced) NRZ data stream, the conditional 
probability of error has the following form: 

*«>.> = Ts e ' k {/S“ s(e ' ) } + 5 ertc {\/S (1 -« + 5"* {\/S (1 - 2 « c “"* ) } 

(37) 

and the variance of the tracking phase error a 2 is given by Eq. (26), where po is defined in Eq. (18) with 


oo 

a = f \H(2TTf)\ 2 S cont (f) df 

— CO 


(38) 


and 


^ tan 2 (m T ) 

Figure 5 shows the curves for SER versus SSNR when data asymmetry is present. 


(39) 


D. Band-Limited Channel 

An additional impairment that contributes to the degradation of the overall performance of the system 
is the ISI caused by the band-limited channel. Band limiting causes interference between successive pulses 
producing the ISI effect, which behaves like an additional random noise. 

If p(t) denotes the pulse shape of the data, and h f (t) denotes the impulse response of the equivalent 
low-pass filter of the RF band-pass filter with bandwidth then the received data can be expressed as 


39 




Fig. 12. Normalized power spectrum for various values 
of data asymmetry. 


oo 

d(t) = dkgit + kT s ) 

k= — oo 

where = ±1 with p = q = 0.5, and g(t) is given by 

g(t) = p(t) * 


(40) 


(41) 


where * denotes convolution. 

The impulse response of an ideal channel h'(t) is given by the inverse Fourier transform of the transfer 
function H'(f ): 


H' 



1 

0 


-B < f < B 
otherwise 


(42) 


resulting in 


ti(t) = 2B 


sin (27r Bt) 
27rBt 


For an ideal filter and a perfect data stream, g(t 4- kT) can be found to be [1] 


( 43 ) 
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(44) 


g(t 4- kT s ) — — 

7T 


Si 


|27rB (t + T s (k + 0 ) | - Si |2tt 5 (t + T s (k - ) } 


where 


Si(x) = J 

o 


sin(u) , 

— du 

u 


(45) 


Figure 13 shows a plot of g(t) versus the normalized time t/T s . Note that the shape of the output is 
dependent on the time- bandwidth product BT S . For BT S » 1, the degradation due to band limiting 
becomes negligible. As BT S approaches 1, the rise and fall times of the output are significant when 
compared to the input, and the output signal is further spread in time. 

To calculate the average probability of error, P e (0 e ) and cr 2 need to be determined. Calculating P e (0 e ) 
exactly is very difficult because one has to take into account all possible combinations of the digits 
d k — ±1 i < |/c| < oo. It is assumed that only a finite number of M pulses before and after d 0 , d 0 = 0, 
are taken into account. That is, only the ISI effects of the M preceding and M subsequent bits are 
considered on the bit under detection. For BT S > 1, the value of 1 < M < 2 is sufficient. The conditional 
error probability may be determined using [1] 


Pe(Oe) 


1 1 
2 2 2M 


£ 


rnrnbmationj) 


erfc 



4- dk A/c 

k^o 


cos (0 e ) 


where the SSNR for this case is given by 


a _ , dt 

N 0 N 0 J 
0 


(46) 


(47) 


and 




T 

f 9{t)9(t 

JO 


+ kT s ) dt 


[ U \9(t )\ 2 
Jo 


dt 


(48) 


The variance of the carrier tracking phase error is given by Eq. (17). Therefore, for 1 < M < 2, the 
average error probability can be obtained by substituting Eqs. (6) and (46) into Eq. (5) and performing 
the numerical integration. The results are shown in Fig. 6. 


III. Combined Effects on PCM/PM/NRZ Receivers 

The practical PCM/PM/NRZ receivers operate in the presence of both an imperfect data stream and a 
band-limited channel. This part of the article studies the combined effects of an unbalanced data stream, 
data asymmetry, and ISI on the SER. The total SSNR degradation of the receivers due to these three 
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Fig. 13. Output response of the ideal filter to NRZ 
pulse for various values of BT 


undesirable effects is not the algebraic sum of the SSNR degradation due to each separate effect found 
in Section II. Therefore, it is necessary to study the combined effects of these three sources on the error 
probability performance. 

For an unbalanced and asymmetric data stream, the dc, continuous, and harmonics-PSD components 
are given by [2] 


Sdc(f) = [ 2 p - (1 - 2 tpt)] 2 6(f) 


(49) 


Scont(f) — T a 


sin 2 ( 71-/7), ) 


+ T, 


[ (*fT a ) 2 

sin 2 (tt/Tj) 


L (*/r ,) 2 


M?t) + «2 (P,Pt,0]+T 3 
[®4 (P,Pt,0 -0-5 (P,Pt)] 


sin 2 (7r fT„0 


{*fT s y 


N(pt,0] 


(50) 


2 00 -1 

ShU) = Y. - mR , ) 

m = l 

respectively, where 


(51) 


®i(Pt) = Pt(l - Pt)[l + 2(1 - p t )) - p 2 


(52) 


42 




MP.Pt.O = {3 Pt +Pt( 1 -Pt)[l + 2(1 — 2p)] } cos 2 (p/T^) (53) 

a 3 (Pt,£) =Pt(l +P ( 2 -P() cos 2 (Trf T a ) +p 2 cos(2tt/T s ^) (54) 

a 4(Pi Pt, 0 =Pt(l -Pt)(l — 2p)(0.5 cos (27 r/T 3 £) - p sin (2 tt/T 3 £)] (55) 

a s(p>Pt) = 0.5p t (l — p t )(l — 2p) (56) 

C(m,p,£) — sin 2 (mn£) [cos 2 (to7t£) — (1 — 2p) 2 sin 2 (m7r£)] (57) 


Note that when p = p t = 1/2 and £ = 0, that is, a perfect data stream, Eqs. (49) through (51) all reduce 
to the result for a perfect NRZ random data stream, Eq. (13). 

Once again, the presence of the two noncoherent tones (the dc component due to the imperfect data 
stream and the carrier tone), both at / c , causes the mean of the phase difference (0q ~ 0) to deviate away 
from zero. The expression for the mean of this phase difference was derived to be 

6 e (mean) = — tan -1 {(tan mr)[(2p — 1) + 2£p(l - p)]} (58) 


Note that this equation reduces to Eq. (25) and Eq. (36) by setting £ = 0 and p — 0, respectively. 

The same approach used in Section II will be used here to determine the average SER. Therefore, 
the simulations will again have to compensate for the phase difference (0o — 0)- The average probability 
of error is given by Eq. (5), and therefore, the expressions for P e {0 e ) and P{0 e ) must be determined to 
evaluate P e . The conditional error probability in the presence of an imperfect data stream and band- 
limited channel is given by [2] 


P e (e e ) =pPr|z(T s ) < ^d 0 = +lj +< 7 Pr|z(T s ) > j,d 0 = -lj (59) 

where the overbar denotes statistical averaging over the joint distribution of the double infinite data 
sequence and the test statistic Z(T S ) is given by 


Z(T S ) = E s 


fc = oo 

±1 + ^2 dk^k(i) 

k = — o c 
k^O 


cos (8 e ) + n(T„) 


(60) 


where ±1 corresponds to do = ±1. It is assumed that the corrupting noise process n(T s ) is a zero-mean 
Gaussian random variable with a variance NqT s /2. The parameter Afc(i) is defined as 


-Mr) 


f 9(t)9i{t + kT s ) dt 
lo 

f W)\ 2 dt 

Jo 


, i = 1,2, 3, 4 


(61) 
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where g(t) is the output of the ideal filter for a perfect data stream given by Eq. (44), and g t (t) for 
i — 1, 2, 3, 4 is defined as [2] 

9i(t + kT s ) = i Si | 2t tB + T s (k + ^) J - Si |2t tB (t + T s (k - 1 - ^)) J (62) 

Si | 2t rfl (t + T s (k + | - Si | 2t xB (« + T s (fc - i + ^ J (63) 

g 3 (t + kT a ) = i 5i|2TB^ + T s ^+^J-5i|27rB^ + T 3 ^-i))| ( 64 ) 

ff4 (« + fcr 3 ) = -ga^ + fcT,) (65) 

The variance of the carrier tracking phase error cr 2 can be obtained using Eq. (26) where 


a= J \H(2nf)\ 2 S cont (f) df (66) 


and 

^ = [2p - (1 - 2 £p t )] 2 tan 2 (m T ) (67) 

Again, a is the interference due to the continuous spectrum component, and I/C is the interference 
caused by the dc component-to-carrier power ratio. The harmonic components caused by asymmetry do 
not interfere with the carrier tracking because of the assumption that 2 Bl « R s - 

The average probability of error can be found by substituting Eqs. (6) and (59) into Eq. (5) and 
performing the numerical integration. The results are shown in Figs. 7 through 9. 


IV. Description of PCM/PM Receiver Blocks 

Figure 14 shows the block diagram of a PCM/ PM receiver. This receiver consists of the test signal 
generator (TSG), the advanced receiver (ARX), and the error counter. The TSG, shown in Fig. 15, 
generates the deep-space spacecraft signal at an intermediate frequency (IF). The TSG’s random data 
block controls the parameter p, and depending on this value, a balanced or unbalanced data stream is 
generated. The data asymmetry block controls the parameter f, producing an asymmetric data stream. 
The Appendix gives a brief description of this block. Setting p = 0.5 and £ = 0 will produce a perfect 
purely random data stream. Setting p / 0.5, f ^ 0, or the combination will result in an unbalanced data 
stream, asymmetric data stream, or a data stream with the combined imperfections, respectively. 
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INPUT PARAMETERS 

SAMPLING RATE, f g : 

CARRIER FREQUENCY, f c 
INITIAL CARRIER PHASE, 8 Q : 
SYMBOL RATE, R g : 

MODULATION INDEX, m f. 

TOTAL POWER/NOISE RATIO, P/N Q : 


5 x 10 5 Hz 
1 xIO 5 Hz 
0.0 deg 
1 xIO 4 Hz 

71.62 deg 
47.455 dB-Hz 


CARRIER PAR AMETERS 

NUMERICALLY CONTROLLED 
OSCILLATOR (NCO) FREQUENCY: 

INITIAL NCO PHASE: 

CARRIER UPDATE RATE: 

ONE-SIDED LOOP BANDWIDTH, B L , CARRIER: 


1 xIO 5 Hz 
0.0 deg 
5x10 5 Hz 
5.0 Hz 



Fig. 14. PCM/PM/NRZ system block diagram as implemented in SPW. 


Othe, TSG parameters include the following (the values shown are the ones used in simulations): 

500 x 10 3 Hz = sampling rate, f 3 
10 x 10 3 Hz = symbol rate, R 3 
100 x 10 3 Hz = carrier frequency, f c 

Odeg = initial carrier phase, 6 0 

71.62 deg = modulation index, m,T (corresponding to 1.25 rad) 
and the total power-to-noise ratio P/N 0 is calculated using 


P 

ivo 


E± 

N 0 


- 101og 10 (sin 2 m T ) + 10 log I0 


R s 


( 68 ) 


where E 3 /N 0 is the SSNR in dB. 

The ARX, shown in Fig. 16, consists of the following blocks: 
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TSG BLOCK PARAMETERS 

SAMPLING RATE, f? 

CARRIER FREQUENCY, 

INITIAL CARRIER PHASE, e Q : 

SYMBOL RATE, R? 

MODULATION INDEX, m f. 


5 x10 s Hz 
1 x 10 s Hz 
0.0 deg 
1 x 10 4 Hz 
71 .62 deg 


TOTAL POWER/NOISE RATIO, PINq. 47.455 dB-Hz 
PROBABILITY OF ZERO, p: 0 5 

NO. OF SAMPLES/SYMBOL 50.0 

DATA ASYMMETRY, ££ 



X 


Fig. 15. Test signal generator (TSG) block diagram. 


(1) The carrier PLL block estimates the incoming carrier phase and frequency and mixes it 
with the input signal. 

(2) The phase imbalance block adds (or subtracts) a phase to its input according to the 
( value tfthe Xse imbalance parameter. The input to this block in simulations is given 

hv Fa (4V therefore depending on the kind of imperfect data stream present, this 
parameter is ^t tode (mean) , as'given by Eqs. (25), (36), or (58), so that by adding this 
phase to the incoming phase 6, the output of the block will have a zero-mean phase error. 
The phase imbalance parameter is set to zero when no phase compensation ,s require , 
that is, when no unbalanced and/or asymmetric data streams are presen . 

(3) The Butterworth low-pass filter controls the presence of the balld "' 1 ^ lt ‘ n ® 

W ting the filter bandwidth B to a value that depends on the product BT S . If no 
limiting is present, the filter bandwidth B is set to 100 kHz. 
u , The ideal clock generates the timing for the sum-dump-hold symbol block The use of 
( > In ideal clock to produce the timing instead of the digital data » 

block was for the purpose of matching the assumption made in theory, and thereto , 

eliminating the loss due to symbol synchronization. 

(5) The sum-dump-hold block outputs the soft symbols. 

Finally, the error counter block compares the soft symbol, of the ARX to the transmitted symbols and 
outputs the number of errors N. 
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CARRIER PARAMETERS 


INPUT PARAMETERS 


NUMERICALLY CONTROLLED 
OSCILLATOR (NCO) FREQUENCY: 

INITIAL NCO PHASE: 

CARRIER UPDATE RATE: 

ONE-SIDED LOOP BANDWIDTH, Bj_, CARRIER: 


1 xIO 5 Hz 
0.0 deg 

5 x 10 5 Hz 
5.0 Hz 


SAMPLING RATE, ^ 5 x 1 0 5 Hz 

CARRIER FREQUENCY, f c 1 x 10 5 Hz 

INITIAL CARRIER PHASE, 0 0 : 0.0 deg 

SYMBOL RATE, Hy 1 x 1 0 4 Hz 

MODULATION INDEX, m f. 71 .62 deg 


TOTAL POWER/NOISE RATIO, P/N Q : 47.455 dB-Hz 



Fig. 16. Advanced receiver block diagram. 


V. Discussion and Simulation Results 


By substituting the expressions for the conditional probability of error P e (^e) 1 “d the 
density function for the phase error P«U into Eq. (5), the SER as a function of SSNR was piottedm 
II 21 for each of the cases discussed above. Using typical operating conditions of m T - 1.25 rad and 
■2B in s = 0 001, these theoretical plots are shown in Figs. 3 through 9 as the continuous curves. The 
computer simulation results are shown as the triangular, circular, and square points for vanables shown 

therein. 


Using the SPW, simulations were performed at 7-, 8-, 9- and 10-dB SSNR (E s /N 0 ) and the corr 
sponding P/N 0 was calculated. The result of each simulation was the number of errors (pro uce y 
the erJ counter as a result of comparing the soft symbols to the transmitted ones). The average error 

probability P e was then calculated using 


P e = 


N_ 

number of iterations/(/,,/ft 3 ) 


(69) 


where / is the sampling frequency in Hz and the fraction (/./*,) is the number of samples per symbol 
The number of iterations must be chosen large enough so that the simulation results have sufficient 

statistics. That is, 


number of iterations = 


100 

SER 


(S 


(70) 
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where SER is the symbol error rate as given by the theory. Finally, P e 
results were compared to the theoretical curves presented in [1,2], ' 


was plotted versus SSNR and the 


A. Unbalanced Data 


To verify the performance of the receiver in the 
were performed for p = 0.5, 0.45, and 0.4. 


presence of an unbalanced data stream, simulations 


As mentioned in Section II B, when p ± 0.5, the phase of the tone caused by the unbalanced data is 
oncoherent with the earner phase, which results in a nonzero-mean phase error 8 e . In order to overcome 

L Pr . e ! erFOr W3S calculated usin S Eq- (25), checked by simulations, and then modified 

so that the resultant phase error is of zero mean. 


Flg wt 3 TA t L heoretlcal and simulation results when no phase modification is made to the phase 
error. When p ^ 0.5 the simulations are in disagreement with theory, and the SSNR degradation in some 
cases exceeds 1 dB. On the other hand, when the phase error is modified, the theoretical and simulation 
results are in good agreement. These results are presented in Table 1 and Fig. 4. It is obvious that as 
P deviates from 0.5 the performance of PCM/PM/NRZ degrades significantly, and that the degradation 
becomes unacceptable when p < 0.45. This is due to the presence of a strong dc component caused by 
the unbalanced data stream at the carrier frequency. The higher the deviation from 0.5, the stronger the 
dc component, and as a result, the worse the degradation. 


Table 1. Simulation data and results for unbalanced data, separate effects. 8 


Probability 
of mark 

E,/N o, 
dB 

. P/N o, 
dB 

No. of 
iterations 
in millions 

No. of 
errors 

Pe 

P e theory 

0.5 

7 

47.455 

6 

105 

8.75 x 10 -4 

7.727 x 10 -4 

0.5 

8 

48.455 

28.2 

126 

2.23 x 10" 4 

1.909 x 10“ 4 

0.5 

9 

49.455 

150 

116 

3.87 x 10“ s 

3.363 x 10~ 5 

0.5 

10 

50.455 

1300 

132 

5.08 x 10“ 6 

3.872 x 10“ 6 

0.45 

7 

47.455 

6.01 

115 

9.57 x 10~ 4 

1.100 x 10“ 3 

0.45 

8 

48.455 

25.2 

129 

2.56 x 10- 4 

3.100 x 10~ 4 

0.45 

9 

49.455 

80.2 

112 

6.98 x 10“ 5 

6.600 x 10- 5 

0.45 

10 

50.455 

500.2 

134 

1.34 x 10~ 5 

1.100 x 10“ 5 

0.4 

7 

47.455 

7 

764 

5.79 x KT 3 

5.400 x 10~ 3 

0.4 

8 

48.455 

10.0 

704 

3.52 x 10- 3 

3.250 x 10 -3 

0.4 

9 

49.455 

15 

612 

2.04 x 10- 3 

2.000 x 10 -3 

0.4 

10 

50.455 

25 

686 

1.37 x 10- 3 

1.500 x 10 -3 

a m = 1.25 rad, R s = 

1 x 10 4 Hz, J,=5x 10 s Hz, B l = 5 Hz, 2 B L /R S = 0.001 


B. Data Asymmetry 

Since the power of the dc component generated by the asymmetric data at f c is much less than the 
power of the carrier tone, the mean of the phase error will be small. The mean was calculated (Eq (36)) 
and measured to be between -1.7 and -5.2 deg for £ between 2 and 6 percent, respectively, which are the 
minimum and maximum values for £ used in the simulations. The degradation due to this nonzero-mean 
p ase error is negligible and, hence, no compensation was done to the phase error. Simulations were 
performed for £ = 2, 4, and 6 percent. The results are shown in Table 2 and Fig. 5. Again the simulation 
results are in good agreement with the theoretical results (within 0.2 dB). The numerical results show 
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that, for data asymmetry less than or equal to 2 percent, the SSNR degradation is on the order of 0.1 dB 
or less, and that this degradation is between 0.2 dB and 0.25 dB for data asymmetry of 6 percent and 
for 1(T 7 < SER < 10" 5 . 


Table 2. Simulation data and results for data asymmetry, separate effects. 8 


Data 

asymmetry, 

percent 

Es/No, 

dB 

Pt/N, 0 , 

dB 

No. of 
iterations 
in millions 

No. of 
errors 

Pe 

P e theory 

2 

7 

47.455 

6.6 

124 

9.39 x 10“ 4 

8.75 x 10" 4 

2 

8 

48.455 

23 

119 

2.59 x 10" 4 

2.20 x 10 -4 

2 

9 

49.455 

130 

139 

5.35 x 10" 5 

3.85 x lO -5 

2 

10 

50.455 

2100 

203 

4.83 x 10“ 6 

4.80 x 10~ 6 

4 

7 

47.455 

6.6 

144 

1.09 x 10~ 3 

9.60 x 10" 4 

4 

8 

48.455 

20 

137 

3.43 x 10 -4 

2.60 x lO' 4 

4 

9 

49.455 

115 

136 

5.91 x 10-® 

4.40 x lO" 5 

4 

10 

50.455 

1900 

253 

6.66 x 10- 6 

5.50 x 10“ 6 

6 

7 

47.455 

6.6 

146 

1.11 x 1(T 3 

1.30 x 10“ 3 

6 

8 

48.455 

18 

125 

3.47 x KT 4 

2.85 x lO" 4 

6 

9 

49.455 

108 

132 

6.11 x l(r s 

4.70 x 10“ 5 

6 

10 

50.455 

800 

138 

8.63 x 10“ 6 

6.30 x lO" 6 

a p = 0.5, m = 1.25 

rad, R s 

II 

X 

© 

Hz, f s = 

5 x 10 5 Hz, 

B L = 5 Hz, 

2 B l /R s = 

0.001. 







C. Band-Limited Channel 

In order to test the effect of the band-limited channel on the overall performance of the system, 
simulations were performed for different values of the time-bandwidth product BT S = 1,2, and 3. As 
expected, the higher the value of the product BT S , the better the performance of the system. The 
simulation results are shown in Table 3 and Fig. 6. The numerical results show that for 10~ 7 < SER < 
10 -5 , the SSNR degradation is in the range of 1 to 1.2 dB for BT S = 1, and less than 0.3 for BT S = 2. 
The theoretical and simulation results are in good agreement. However, the simulations are a little worse 
than the theoretical results. This is because the theoretical results were obtained for the case when the 
ISI is caused by two adjacent pulses, that is, two pulses before and two pulses after the current pulse is 
considered in the SER calculation. 

D. Combined Effects 

To test the behavior of PCM/PM/NRZ receivers in the presence of the combination of the three 
undesirable effects, simulations were performed for different values of p, £, and BT S . One of the parameters 
was varied as the other two remained constant. Since data imbalance and asymmetry were always present, 
all simulations required compensation for the phase error 6 e (mean) (Eq. (58)) so that the result is a zero- 
mean phase error. 

Figure 7 plots the SER as a function of SSNR for a fixed data asymmetry f of 2 percent and BT S = 3 
with p, probability of mark, as a parameter. The simulation results are also shown in Table 4, and are in 
good agreement with the theory. The results indicate that, for mr = 1-25 rad and 2 Bi/R$ = 0.001, the 
SER degrades seriously as p deviates from 0.45. 
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Table 3. Simulation data and results for a band-limited channel, 
separate effects* 


BT a 

E./No, 

dB 

P/No, 

dB 

No. of 
iterations 
in millions 

No. of 
errors 

P C 

P e theory 

1 

7 

47.455 

6 

226 

1.88 x 10- 3 

1.80 x 10~ 3 

1 

8 

48.455 

11 

141 

6.41 x lO -4 

5.80 x 10- 4 

1 

9 

49.455 

30 

94 

1.57 x 10- 4 

1.70 x 10- 4 

1 

10 

50.455 

160 

88 

2.75 x 10“ 5 

3.30 x 10“ 5 

2 

7 

47.455 

6 

150 

1.25 x 10“ 3 

9.20 x 10“ 4 

2 

8 

48.455 

21 

146 

3.48 x 10" 4 

2.50 x 10~ 4 

2 

9 

49.455 

105 

143 

6.81 x 10- 5 

4.83 x 10 -5 

2 

10 

50.455 

800 

137 

8.56 x 10“ 6 

6.50 x 10" 6 

3 

7 

47.455 

6 

125 

1.04 x 10 -3 

8.60 x 10~ 4 

3 

8 

48.455 

22.6 

131 

2.90 x 10~ 4 

2.30 x 10“ 4 

3 

9 

49.455 

114 

126 

5.53 x 10" 5 

4.40 x 10~ 5 

3 

10 

50.455 

800 

124 

7.75 x 10“ 6 

5.60 x 10~ 6 

R 771 = 

Bl 

1.25 rad, probability of mark 
= 5 Hz, 2 B l /R s =0.001. 

= 0.5, fi, = lx 10 4 

Hz, /, = 5 x 10 s Hz, 


Table 4. Simulation data and results for various probabilities of mark, 
combined effects. 8 


Probability 
of mark 

E s /N o, 
dB 

P/N o, 
dB 

No. of 
iterations 
in millions 

No. of 
errors 

Pe 

P € theory 
(approximate) 

0.45 

7 

47.455 

7.2 

149 

1.03 x 10~ 3 

7.80 x 10“ 4 

0.45 

8 

48.455 

25.2 

137 

2.72 x 10- 4 

2.00 X 10“ 4 

0.45 

9 

49.455 

80.2 

78 

4.86 x 10“ 5 

3.60 x 10- 5 

0.45 

10 

50.455 

500.2 

57 

5.70 x 10“ 6 

5.40 x 10~ 6 

0.4 

7 

47.455 

7.2 

467 

3.24 x 10“ 3 

2.20 x 10~ 3 

0.4 

8 

48.455 

25.2 

711 

1.41 x 10~ 3 

1.25 x 10~ 3 

0.4 

9 

49.455 

80.2 

983 

6.13 x 10“ 4 

6.50 x 10- 4 

0.4 

10 

50.455 

500.2 

4.23 x 10 3 

4.23 x 10- 4 

3.85 x 10~ 4 

0.35 

7 

47.455 

7.2 

2.80 x 10 3 

1.94 x 10“ 2 

1.50 x 10“ 2 

0.35 

8 

48.455 

25.2 

6.41 x 10 3 

1.27 x 10 -2 

1.25 x 10~ 2 

0.35 

9 

49.455 

80.2 

1.89 x 10 4 

1.18 x 10~ 2 

1.00 x 10- 2 

0.35 

10 

50.455 

500.2 

8.72 x 10 4 

8.72 x 10“ 3 

9.00 x 10- 3 

“Data asymmetry = 2 percent, BT 3 = 3, m = 1.25 rad, R s = 1 x 10 4 Hz, 
f, = 5x 10 5 Hz, B l = 5 Hz, 2 B L /R S = 0.001. 


Table 5 shows the simulation results obtained for various values of data asymmetry £ with BT S 3 
and p = 0.45. As shown in Fig. 8, the simulations are in good agreement with the theory. It is also 
obvious that PCM/PM/NRZ is not sensitive to data asymmetry since the SSNR degradation^ between 
0.1 and 0.5 dB when £ varies between 2 and 6 percent and the SER is between 10 and 


Table 5. Simulation data and results for various values of data 
asymmetry, combined effects. 


Data 

asymmetry, 

percent 

E s /N o, 
dB 

No. of 

°* iterations 
dB .... 

in millions 

No. of 
errors 

Pe 


P e theory 
(approximate) 

2 

7 

47.455 

7.2 

149 

1.03 x 

10~ 3 

7.80 x 

icr 4 

2 

8 

48.455 

25.2 

137 

2.72 x 

IO" 4 

2.00 x 

10“ 4 

2 

9 

49.455 

80.2 

78 

4.86 x 

10~ 5 

3.60 x 

10 -5 

2 

10 

50.455 

500.2 

57 

5.70 x 

10~ 6 

5.40 x 

10 -6 

4 

7 

47.455 

7.2 

160 

1.11 X 

10“ 3 

9.70 x 

io- 4 

4 

8 

48.455 

25.2 

155 

3.08 x 

10- 4 

2.70 x 

ltr 4 

4 

9 

49.455 

80.2 

95 

5.92 x 

10~ 5 

4.80 x 

10~ 5 

4 

10 

50.455 

500.2 

78 

7.80 x 

10~ 6 

6.50 x 

icr 6 

6 

7 

47.455 

7.2 

168 

1.17 x 

10" 3 

1.10 x 

io- 3 

6 

8 

48.455 

25.2 

168 

3.33 x 

10- 4 

2.95 x 

io -4 

6 

9 

49.455 

80.2 

98 

6.11 x 

ur 5 

5.50 x 

io - 5 

6 

10 

50.455 

500.2 

82 

8.20 x 

10 -6 

7.75 X 

10" 6 

a Probability of mark = 0.45, BT S — 3, 
f s = 5 x 10 5 Hz, B l = 5 Hz, 2B l /R s = 

m = 1.25 rad, 
0.001. 

R s - 

: 1 X 10 4 HZ, 


Table 6 and Fig. 9 illustrate the SER performance in the presence of a band- limiting channel for 
p = 0 45 and £ = 2 percent with BT S as a parameter. As shown, the simulations are in good agreemen 
with the theory, and for BT S = 3, the SSNR degradation is on the order of 0.4 dB or less when the SER 

is between 10" 4 and 10“ 7 . 

The numerical results prove that the total SSNR degradation due to the three undesirable effects is 
not the algebraic sum of the SSNR degradation due to each separate effect As an example when th 
SER is 1CT 5 , the SSNR degradation when p = 0.45, £ - 2 percent, and BT S 3 ( lg. ) ' . . ’ 

whereas, the algebraic sum of the SSNR degradations due to each separate effect (Figs. 4 through 6) 

about 0.6 dB. 


VI. Conclusion 

This article studied, by computer simulations, the separate and combined effects of unbalanced data, 
dautymmtrj! and 1 Unlimited channel on ,he performance of a PCM/PM/NRZ r»«M the 
simulation results were in good agreement with the theoretical results presented in [ , ]• ’ 

mathematical models presented in [1,2] can be used to predict the performance 0 !' 

receivers PCM/PM/NRZ was shown to be most sensitive to the imbalance between +1 s and i s in 
the data stream, as the performance degradation became unacceptable when p < R45 and 1 least ^sensi i ive 
to data asymmetry. For BT S = 3, the SER performance was shown to be acceptable for both near- Earth 

and deep-space missions. 
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Table 6. Simulation data and results for various values of BT 8 , 
combined effects.® 


BT S 

E./No, 

dB 

P/N o, 
dB 

No. of 
iterations 
in millions 

No. of 
errors 

Pe 

P e theory 
(approximate) 

3 

7 

47.455 

7.2 

149 

1.03 x 10“ 3 

7.80 x 10 -4 

3 

8 

48.455 

25.2 

137 

2.72 x 10" 4 

2.00 x 10“ 4 

3 

9 

49.455 

80.2 

78 

4.86 x 10“ 5 

3.60 x 10~ 5 

3 

10 

50.455 

500.2 

57 

5.70 x 10‘ 6 

5.40 x 10 -6 

2 

7 

47.455 

7.2 

168 

1.17 x 10~ 3 

8.50 x 10" 4 

2 

8 

48.455 

25.2 

158 

3.13 x 10“ 4 

2.20 x 10" 4 

2 

9 

49.455 

80.2 

92 

5.74 x MT 5 

3.90 x 10" 5 

2 

10 

50.455 

500.2 

70 

7.00 x 10~ 6 

6.60 x 10~ 6 

1 

7 

47.455 

7.2 

222 

1.54 x 10“ 3 

1.25 x 10" 3 

1 

8 

48.455 

25.2 

244 

4.84 x 10 -4 

3.75 x 10~ 4 

1 

9 

49.455 

80.2 

177 

1.10 x 10“ 4 

1.20 x 10" 4 

1 

10 

50.455 

500.2 

207 

2.07 x ur 5 

2.60 x 10" 5 

a Probability of mark = 0.45, data asymmetry = 2 percent, m 
1 x 10 4 Hz, /, = 5 x 10 6 Hz, B l = 5 Hz, 2 B L /R S = 0.001. 

= 1.25 rad, R s = 


Another modulation scheme that is of interest to CCSDS is PCM/PM/Bi-0, which is also known 
to be one of the most efficient modulation schemes in terms of bandwidth occupancy as compared 
to PCM/PSK/PM [4], Mathematical models have been developed to predict the performance of 
PCM/PM/Bi-0 receivers [1,2], and these models are currently being verified by members of CCSDS 
and the results will be reported later. 
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Appendix 

Data Asymmetry Block 


The data asymmetry block outputs NRZ asym mykm ^LlTand 

ZZ - - ° f - — " 8i " 6 

dk - dfc-i (A-l) 


trans = 


iii A H is the previous symbol value, and therefore, this yields 
where d k is the present symbol value and d k i +1 transition occurs, and a 0 when no 

a -1 when a +1 to -1 transition occurs, a +1 when a 1 to ^ Q if frans = +1 or 0, and € 

transition occurs. The block then determmes a re samples per symbol for the input x, P 

if *»»• - -1. ^ere € deno^s daU = et T f ther^ ^ ^ then , for i « 0, . < N; 

denotes the past sample value, C the current a = 

if i < Ti, then y = P; otherwise, if Ti < r < N, the J/ 
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Recent improvements in computational capability and DSN t^h^i u 
newed interest in examining the dos sihilitv nf , - Y d ° bN 1 “ nolo Sy have re- 

to navigate interplanetary spacecraft ThJ S1D S one-way Doppler data alone 
stondjd dideJcedcoTni Zplfor “ <*• 

A covariance analysis, which « T^lTf “ T 

one-way Doppler data is DerfnrmL an w ac J obtail mble by combinations of 

single-station data, the position dispersions are roughly TordZZm7Zh77del 

re 7 **, ~ to 


I. Introduction 

With increasing emphasis on controlling the costs of rWn ™ a ,. • • 

examined that decrease the costs of the spacecraft itself n P P v m,ssl °ns, several options are being 
coherent mode; that is, the spacecraft^ TeTStar^a i.™! ”7 ° Pt '°” jf 10 ^ in » »°» 

carrier signal. Historically, one-way Doooier data h a ^ponder capable of coherently returning a 
instability of spaceborne oscillators, the useof S-band (2 3 < GHz^fr USe ^ s ? le data ^Po due to the 
sources that could not be adequately modeled, Howevefthh the Xni‘1’K“n c “ res e‘>" di »8 error 
more sophisticated modeling ability the possibility nf " ^ v ® nt of h 'gh-speed workstations and 

article assesses the navigation pirn ' 77 “ L 6 “ iy D ° Ppler ls bein * cammed. This 
e^y missions. As a reUnEXZ^ f ” « "-U- 

trajectory were used to perform the analysis r„i r ’ th e Mars Pathfinder spacecraft model and 
data formulated as standard differenced-count Dop^ler^wh^h bet f W6en results em P lo y in g Doppler 
as accumulated carrier phase (which yields a distance L ** ^ ^ & frequency measurement) as well 
Combinations of one-way data obtained simultaneously at two^-ff 6 " ’ USUal,y glven in terms of cycles), 
produce an angular type measurement) and single-statkm one-wa hT* Stat “ nS and then differenced (to 
may satisfy future mission requirements. y a a are shown to produce results that 
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